Integrated bioinformatics and network pharmacology to identify the therapeutic target and molecular mechanisms of Huangqin decoction on ulcerative Colitis

Huangqin decoction (HQD) is a Traditional Chinese Medicine formula for ulcerative colitis. However, the pharmacology and molecular mechanism of HQD on ulcerative colitis is still unclear. Combined microarray analysis, network pharmacology, and molecular docking for revealing the therapeutic targets and molecular mechanism of HQD against ulcerative colitis. TCMSP, DrugBank, Swiss Target Prediction were utilized to search the active components and effective targets of HQD. Ulcerative colitis effective targets were obtained by microarray data from the GEO database (GSE107499). Co-targets between HQD and ulcerative colitis are obtained by Draw Venn Diagram. PPI (Protein–protein interaction) network was constructed by the STRING database. To obtain the core target, topological analysis is exploited by Cytoscape 3.7.2. GO and KEGG enrichment pathway analysis was performed to Metascape platform, and molecular docking through Autodock Vina 1.1.2 finished. 161 active components with 486 effective targets of HQD were screened. 1542 ulcerative colitis effective targets were obtained with |Log2FC|> 1 and adjusted P-value < 0.05. The Venn analysis was contained 79 co-targets. Enrichment analysis showed that HQD played a role in TNF signaling pathway, IL-17 signaling pathway, Th17 cell differentiation, etc. IL6, TNF, IL1B, PTGS2, ESR1, and PPARG with the highest degree from PPI network were successfully docked with 19 core components of HQD, respectively. According to ZINC15 database, quercetin (ZINC4175638), baicalein (ZINC3871633), and wogonin (ZINC899093) recognized as key compounds of HQD on ulcerative colitis. PTGS2, ESR1, and PPARG are potential therapeutic targets of HQD. HQD can act on multiple targets through multi-pathway, to carry out its therapeutic role in ulcerative colitis.

Estrogen receptor PPARG Peroxisome proliferator-activated receptor gamma Ulcerative Colitis (UC) is one subtype of inflammatory bowel disease (IBD) of unknown aetiology that affects the colorectum, often leads to bloody diarrhoea and abdominal pain 1 . Multiple pathogenesis, such as environmental factors, genetically susceptibility and gut microbiome dysfunction, have been implicated 2 . The mucosal inflammation initiating in the rectum and extending proximally in the colon, and limited to the mucosal layer, contrast to Crohn's disease (CD) is involving all layers of the bowel wall 1 . Owing to intestinal homeostasis dysfunction and intestinal barrier disruption, may result in long-term relapsing and remitting mucosal inflammation 3 . Patients with ulcerative colitis experience unpredictable flares-up may present fever, tenesmus, anemia, and possibly intestinal perforation, which is an afflicted condition of an individual 4,5 . Despite 5-ASA, rectal corticosteroids, and biological drugs are the main treatment of ulcerative colitis, approximately 15% of patients may require restorative proctocolectomy with ileal pouch-anal anastomosis (IPAA) for medically refractory disease or to treat colonic neoplasia 6,7 . However, with the incidence and prevalence of ulcerative colitis is increasing worldwide, governments should be informed of the management and economic burden of this disease. Therefore, it is urgent to seek more therapeutic armamentariums with fewer adverse side effects for ulcerative colitis. Traditional Chinese Medicine (TCM) has almost 5000 years of history, which creates many treatments for disease. Huangqin decoction (HQD), a traditional herbal formula from "Shang Han Lun" written by Zhang Zhongjing in Han Dynasty, which is used for "heat-dampness" ZHENG (resulting in fever, abdominal pain, and hematochezia) for over 1800 years. HQD consists of four herbs: Scutellariae Radix (SR, Huangqin), Paeoniae Radix Alba (PRA, Baishao), Jujubae Fructus (JF, Dazao), Glycyrrhiza glabra L. (GL, Gancao). Recently, some fundamental analysis has proved that HQD may regulate the gut microbiota, suppress Ras-PI3K-Akt-HIF-1α and NF-κB pathways 8 , target ESR1 and EGFR to relieve endothelial dysfunction 9 to treat ulcerative colitis. However, pharmacology and molecular mechanisms remain unclear.
In this study, we combined network pharmacology and bioinformatics to reveal the core target and pathways that are participated in the pathogenesis of HQD against ulcerative colitis and screen for molecular targets of HQD for the treatment of ulcerative colitis. The workflow is shown in Fig. 1.

Protein-protein interaction (PPI) network. Standardize gene names from Uniprot Knowledgebase
(UniprotKB, http:// www. unipr ot. org), selecting "Homo sapiens". Draw the Venn diagram (http:// bioin forma tics. psb. ugent. be/ webto ols/ Venn/) for co-targets between HQD and DEGs and shown as a heatmap. Built on the co-targets result, the PPI network was conducted by Search Tool for the Retrieval of Interacting Genes (STRING, https:// www. string-db. org/ 14 , then visualized by Cytoscape3.7.2 15 . The core target was topological analyzed using the CytoNCA plugin, limited core targets over 1 time medians of degree, betweenness, and closeness, respectively. GO and KEGG pathway enrichment analysis. Metascape (http:// metas cape. org/ 16 serves as a gene annotation and analysis resource platform, limited "H.sapiens" and P < 0.01 for GO and KEGG pathway enrichment analysis 17 Table S1). Among 161 active components, 468 effective targets were selected by the Swiss Target Prediction and DrugBank database. Twenty core components as interactions over 100 were listed in Table 1.

Identification of DEGs and co-targets.
According to the GEO2R web tool, a total of 1542 DEGs were identified, which consist of 1024 upregulated and 518 downregulated genes ( Fig. 3 and Supplement Table S2). Venn analysis between 1542 DEGs and 468 HQD targets, had got 79 co-targets (Fig. 4A). Hierarchical cluster analysis of 79 co-targets between ulcerative colitis mucosal samples and normal controls in GSE107499 (Fig. 4B).
PPI network and core target. In the PPI network (      (Table 2 and  Supplement Table S3).
GO and KEGG pathway enrichment analysis. The GO and KEGG pathway enrichment analysis of 79 overlapping targets were performed relying on Metascape platform according to the P < 0.01. There were 818, 26, and 66 GO terms related to biological processes, cell components, and molecular functions, respectively (Fig. 6). Inflammatory response, collagen metabolic process, and leukocyte migration was the significantly enriched for BP. For Membrane raft, extracellular matrix, and side of membrane was the significantly enriched for CC. Oxi-  www.nature.com/scientificreports/ doreductase activity, endopeptidase activity, and protein tyrosine kinase activity was the significantly enriched for MF. In addition, 184 pathways were contained in KEGG pathway enrichment analysis, mainly including TNF signaling pathway, IL-17 signaling pathway, Rheumatoid arthritis, Cytokine-cytokine receptor interaction, Th17 cell differentiation (Fig. 7). Molecular docking. Molecular docking was applied to analyse the binding of the core targets (IL6, TNF, IL1B, PTGS2,ESR1, and PPARG) and all 19 active compounds of HQD (except for MOL001490, 3D conformer generation is disallowed on PubChem since too flexible), respectively. The binding energy score was all smaller than − 5.0 kcal·mol −1 , and the results are listed in Table 3.
Furthermore, according to the ZINC15 database, quercetin (ZINC4175638), baicalein (ZINC3871633), and wogonin (ZINC899093) as the stable compound recognized as key compounds. From Table 3, the binding energy of 5F19 and 1PRG was lower than others, which indicated that PTGS2 and PPARG were the target of HQD on ulcerative colitis.

Discussion
Given its unclear pathogenesis, potentially progressive, and debilitating disease course, the therapeutic goals for ulcerative colitis have changed from treating symptoms to mucosal healing. However, conventional treatments are mainly restricted to control inflammation and clinical symptoms, and hard to meet the expectations of long-term management. Therefore, combination therapy is the urgent strategy for ulcerative colitis 18 . Profound knowledge of genetics, microbiota or cytokine imbalances in the bowel and its function in inflammation will sustain the outcome of possible therapeutic approaches if we act on those targets. In addition, the development of new ligands with therapeutic efficacy on the bowel would be a potential avenue in the search for treatments or prevention of ulcerative colitis. To date, with widely using TCM in different diseases, clinical effectiveness has been demonstrated. As for this concern, we made a combination of network pharmacology and bioinformatics might contribute to a further understanding of the disease aetiology and find novel potential therapeutic way. www.nature.com/scientificreports/ In this study, 20 core components with 486 targets of HQD and 1542 DEGs (GSE107499) of ulcerative colitis were examined. 79 co-targets which indicated that the targets of HQD against ulcerative colitis may be through IL6, TNF, IL1B, PTGS2, ESR1,PPARG, which have a high degree in PPI network. IL6, TNF, IL1B have been proved that have a strong impact on the pathological process, as for docking test, we found PTGS2 and PPARG have the lower energy than others, and we focus on more details in PTGS2, ESR1, and PPARG. ESR1, a nuclear hormone receptor, is involved in the regulation of eukaryotic gene expression and affects cellular proliferation and differentiation in target tissues by the steroid hormones and their receptors. ESR1 has mutual transrepression with NF-κB, e.g. decreases NF-κB DNA-binding activity and inhibits NF-κB-mediated transcription 19 . Also, it can involve activation of NOS3 and endothelial nitric oxide production 20 . Dysregulation of estrogen receptors in the intestinal mucosa of ulcerative colitis indicates that estrogen signaling may play a role www.nature.com/scientificreports/ in the local immune response and maintain epithelial homeostasis in a gender-and age-dependent manner 21 , even though there is no evidence in difference prevalence in male and female. PPARG (PPARγ), a new coming experiment targeting the treatment of ulcerative colitis, is one of the isoforms of nuclear PPAR receptors that belong to the steroid receptor superfamily. It can be bound to DNA-specific PPAR response elements (PPRE) and modulates the transcription of its target genes, so that controls the peroxisomal beta-oxidation pathway of fatty acids. In addition to playing a role in the regulation of cardiovascular circadian rhythms adipocyte differentiation, and glucose homeostasis, PPARγ also acts as a critical regulator of gut homeostasis by suppressing NF-κB-mediated proinflammatory responses [22][23][24] .
Presence of PPARγ can be observed throughout the gastrointestinal tract: mesenteric adipocytes, macrophages, and epithelium, most prevalently in the more differentiated layer of epithelial cell, have the highest levels of PPARγ mRNA than other tissue [25][26][27][28] . However, disease onset of ulcerative colitis provokes inflammation of the colonic intestinal wall, showing significant downregulation of PPARγ, compared with remission mucosal tissue [29][30][31] . Increased transcriptional activity of PPARγ was linked to decreased intestinal inflammation in DSSinduced colitis mice 32 .
Besides, a present study has proved that PPARγ can be activation by COX-2 products (15d-PGJ2-G, a chemical metabolite of PGD2-G) to perform anti-inflammatory in ulcerative colitis 33 . PTGS2 (COX-2) is a key gene with a particular role in the inflammatory response and expresses in epithelial cells and mononuclear cells in IBD 34,35 . Prostaglandins are critical to the process of mucosal healing in the colonic tract 36 .
Moreover, core components of HQD act in a multi-component and multi-target way. Quercetin, acacetin, wogonin, norwogonin, baicalein, salvigenin, oroxylin a were key active components of HQD in the treatment of ulcerative colitis. Quercetin has an antioxidant effect, suppresses proinflammatory mediators' release and expression of inflammatory proteins 37 . Acacetin by inhibiting inflammation and regulating the intestinal microbiota in DSS-induced colitis mice 38 . Wogonin is quickly converted into baicalin and various flavones in metabolism and may act in the PI3K-Akt pathway, AMPK pathways, and inhibition of telomerase activity 39 . Baicalin has been confirmed to inhibit IL-33 and NF-κB p65 expression, increase IκB-α levels, regulate Th17/Treg balance and gut microbiota to act antiinflammation 40,41 . Oroxylin A as the main component extracted from Scutellariae radix can decrease the expression of inflammatory cytokines IL-6 and IL-1β, and attenuated IL-6/STAT3 signaling pathway in ulcerative colitis mice 42 . Norwogonin can increase Bax and decrease in Bcl-2 levels, induction of autophagy along with triggering G2/M phase cell cycle arrest to suppress the colon cancer 43 .
In GO and KEGG pathway enrichment analysis, HQD may through several cellular components, such as membrane raft, extracellular matrix, and receptor complex to active in molecular functions, including oxidoreductase activity, endopeptidase activity, and protein kinase activity to get involved in biological processes, such as TNF signaling pathway, IL-17 signaling pathway, Th17 cell differentiation, and NF-kappa B signaling pathway.
Recover intestinal mucosal homeostasis faces a particularly challenging change, because it is against infection, commensal organisms, and unbalance T cells differentiation 44 . Th17 is the most synthetic T cell, plays a www.nature.com/scientificreports/ vital role in gut between homeostasis and inflammation, and is thought to play a pathogenic role in IBD 45,46 . In a previous study, Th17 cell differentiation can be driven by lnc-ITSN1-2 and correlate with increased IBD risk, activity, and inflammatory cytokines 47 . NF-κB signaling is often dysregulated resulting in overzealous inflammation and distinctly into two pathways, canonical or non-canonical 48 . Different from rapid, acute inflammation is canonical NF-κB signaling pathway, non-canonical NF-κB signaling has become more evident in maintaining immune system homeostasis in the gut, which regulates T cells differentiation and function 49,50 . In sum, HQD mainly acts into multi-component, multi-target, multi-pathway in treating ulcerative colitis.

Limitations
Further pharmacological experiments are needed in order to validate the therapeutic mechanism of HQD on ulcerative colitis.

Conclusion
In this study, we combined bioinformatics analysis, network pharmacology, and molecular docking to reveal the mechanisms of HQD against ulcerative colitis. HQD may be through multi-components and multi-pathway in treating ulcerative colitis.